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Abstract 

Nonmagnetic impurities in d-wave superconductors cause strong local sup- 
pressions of the order parameter. We investigate the observable effects of 
the scatterigng off such suppressions in bulk samples by treating the order 
parameter "hole" as a pointlike off-diagonal scatterer treated within a self- 
consistent t-matrix approximation. Strong scattering potentials lead to a 
finite-energy spectral feature in the d-wave "impurity band", the observable 
effects of which include enhanced low-temperature microwave power absorp- 
tion and a stronger sensitivity of the London penetration depth to disorder 

than found previously in simpler "dirty" d-wave models. 
PACS Numbers: 74.25.Dh, 74.25.Nf, 74.62.Dh 
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Introduction. The current consensus regarding the d-wave symmetry of the order pa- 
rameter in the hole-doped cuprates rests to a considerable extent on the success of simple 
models of the effect of disorder on unconventional superconductors. [|1] Theories based on 
the impurity t-matrix approximation |2|-|J suggest that disorder gives rise to an "impurity 
band" of extended states near the Fermi level which dominates the low-temperature be- 
havior of the system. Although published data on superconducting properties continue to 
exhibit a certain variability such apparent discrepancies can frequently be accounted for 
within this picture by allowing for the variation in sample quality. Systematic substitution 
of Copper with nonmagnetic planar defects such as Zn or Ni leads to similar effects. The 
observed strong temperature dependence of transport properties of nominally pure samples 
with nearly optimal T c gives rise to the further conclusion of near-unit arity limit scattering 
by most simple defects in these systems. 

Several puzzles and discrepancies must be understood before the problem of disorder 
in the superconducting state of the high-temperature superconductors (HTSC's) may be 
considered solved even at a purely phenomenological level. The simplest theory of low- 
temperature transport in a d-wave superconductor |JJ predicts an u> — > conductivity 
a ~ (T 00 + a T 2 , where cr o — ne 2 / (mTrA ) is a "universal" constant dependent on the scale 
set by the order parameter maximum Ao but independent of the relaxation rate 1/r for 
weak disorder. || Experiments on untwinned samples clean enough to exhibit the linear-T 
magnetic penetration depth characteristic of a pure d-wave state down to a few degrees 
are consistent with the universal residual conductivity value, but exhibit a conductivity 
temperature variation 8a ~ T or even sublinear temperature dependence. || 

Further difficulties include the rather slow measured rate at which the critical temper- 
ature itself is supressed by impurities, (for planar defects, roughly a factor of two slower 
than predicted by the simplest theory) as well as a considerably more rapid decrease of 
the low-T superfluid density p s with disorder than predicted. ]7]|| The two effects together 
result in a ratio T c j p s which is considerably larger than predicted by "dirty d-wave" theory 
in many HTSC samples, as pointed out recently by Franz et al. [JToJ] These authors argued 



that the usual theory, which assumes a spatially homogeneous disorder-averaged order pa- 
rameter (OP) Ak, neglects an important physical effect, namely the supression of the true 
order parameter Aj c (R) around the positions of the defects. These static fluctuations are 
expected to be large since the order parameter is d-wave in character, and highly local since 
the order parameter varies (crudely speaking) on a scale of the coherence length £, which 
at low T is of order a few A in HTSC's. In a full numerical solution to the Bogoliubov-de 
Gennes (BdG) equations for a d-wave system with dilute strong potential scatterers, Franz 
et al. indeed observed an enhancement of T c j p s consistent with experiment. Very recently 
ljj , a more analytic approach to the effect of order parameter suppression by Zhitomirsky 
and Walker found a reduction of the T c -suppression in rough agreement with experimental 
observations. Both approaches work best close to T c . On the other hand, the discrepancy 
with experiment in the T-dependence of the transport properties is most evident at low 
temperatures. Clearly, a complementary approach for low T is needed. 

In the interest of constructing a practical theory of transport in the HTSC's including 
these order parameter suppressions, we have attempted to exploit the short range nature of 
the order parameter supression (at low T) by replacing the true, self-consistently determined 
Ak(R) near the impurity site by a pointlike order parameter "hole" which acts as an off- 
diagonal (in the Nambu-Gorkov matrix sense) potential scatterer for electrons. The weight 
of the off-diagonal perturbation is then determined by the bare impurity potential and 
by the solution to the one-impurity problem for the R-integrated OP suppression. This 
approximation makes sense at low temperatures, mainly because the zero T coherence length 
which controls the range of the OP suppressions is at most a few atomic lengths for the 
HTSC's. With this approximation, one can perform disorder averages in the usual way and 
obtain a generalization of the very flexible t-matrix approximation involving a translationally 
invariant, effective medium. 

Having outlined the method, we proceed to show that the failures of the "dirty d-wave" 
model when compared to experiments on low-T microwave conductivity and the dependence 
of the absolute penetration depth on disorder are at least partially cured. The results for the 
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temperature dependence of A and other thermodynamic quantities does change very little. 



A preliminary account of this research has appeared in Ref. |L2 . 

Local order parameter supression by single impurity. Several authors have solved the BdG 
equations or equivalent for the local structure of the d-wave Ak(R) around a nonmagnetic 
impurity. In our approach, we make the usual BCS assumption that the pairing potential 
is separable, Vk,k' = V<&(k)<&(k'), and has d-symmetry, with $d(k) = -\Z2cos20 normalized 
over a model circular Fermi surface. We further neglect so-called "leading loser" components 
of the pair interaction, e.g. subdominant pairing channels which have been shown to lead 
to additional fine structure in the order parameter around the impurity site despite being 
energetically forbidden in the bulk. 

The bare impurity itself is described for simplicity by a ^-function potential in real space, 
?7(R — Rj mp ) = U 5(Tl — Rj mp )r 3 , where the Tj are the Pauli matrices in particle-hole space. 
Initially, one might try to find the spatial variation of the order parameter, which is given 
by the BCS gap equation after subtraction of the bulk limit. The Fourier transform of the 



(static) order parameter fluctuation <5Ak(q) is then determined |T^JT^] by the single- impurity 
t-matrix T(p, p'), 

5A k (q) = V$ d (k)T E <^( k ') >< (!) 

w k' 

Tr { yG (k' + q/2)T(k' + q/2, k' - q/2)G (k' - g/2)} . 

(2) 

where Go is the matrix Green function of the pure system. The t-matrix is given as usual 
by 

T( P , p') = t>(p, p') + £ t>( P , p")G (p")T(p", p '), ( 3 ) 

p" 

where T(p, p') and f/(p, p') are Fourier transforms with respect to the electronic momenta. 
In the usual "dirty d-wave" theory, the t-matrix is taken independent of momentum, T = 
T(u), for the case of isotropic scatterers. Here, we explicitly account for the fact that 
electrons moving in the neighborhood of the impurity feel an effective one-body potential due 



not only to the bare impurity but to the order parameter modification about the impurity. 
The effective impurity potential is therefore 

f>k(q) = f> T- 3 + 5A k (q)r 1 (4) 

where k = (p + p')/2 is the momentum conjugate to the "fast" relative motion of the 
electron pair, while q = p — p' is the momentum conjugate to the slowly varying center of 
mass position. Solving Eqs. self-consistently is equivalent to solving the BdG equations 
for a single impurity, and will yield a momentum-dependent t-matrix in the d-wave case. 
The associated order parameter fluctuation <5Ak(q) has d-wave symmetry in real space about 
the impurity site. [13| 



Our objective instead is to develop a formalism for calculating observables in the presence 
of a finite concentration of impurities. To this end, we neglect the q-dependence of 5Ak(q) 
in Eq. £| equivalent to replacing the full <5Ak(R — Rjmp) with <5Ak(q = 0) <5(R — R; mp ). 
This is certainly a reasonable approximation in the cuprates at temperatures sufficiently far 
below T c , since the very short T = coherence length £ on which the order parameter varies 
becomes comparable to the lattice constant. We then solve Eq. |2| for <5A k (q = 0) under the 
continued assumption of negligible subdominant pair component, such that <5A k (q = 0) = 
Sd cos 20. We point out that no new parameter has been introduced into the theory by this 
procedure, since 5j is driven by the impurity scattering strength Uq. 
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FIG. 1. Order parameter scattering strength vs. temperature. The dot-dashed and dashed 
lines are the first order (FO) result [15| for 5d with Uq = 10 and U = 20, respectively, the symbols 
and solid lines are the self-consistent results for C7o=10,20,50,1000, as indicated. For this example, 



100T C and A D = 3T C . The pairing potential and the cutoff are chosen such that the 



BCS-expression for T c is the unit of energy, T c = 1.14cj exp (—1/N \V\) (u> = 30, \V\ ~ 28.31) 

We begin by solving Eqs. @-|4] numerically under the approximations outlined above. 
In Fig. |l] we show the temperature dependence of the self-consistent evaluation of 3d in 
comparison with the result obtained when the T\ part of the t-matrix is iterated only once. 
j~5| Close to T c , the first order (FO) result diverges like (1 — T /T c )~ l l 2 whereas the self- 



consistent result vanishes with the bulk order parameter Ao like (1 — T /T c ) 1 ^ 2 . For low 
T and the smallest Uq = 10 the first order and the self-consistent result agree well, as 3d 
is sufficiently small. But already for Uq = 20 the discrepancies are obvious even for zero 
temperature. For small Uq (and low T) the first order result underestimates \Sd\, however, 
for large Uq it overestimates \Sd\ (not shown in the figure). Note that Sd can be quite large 
in the unitarity limit Uq — > oo. However, it can never be infinite as its f/o^dependence 
saturates like Uq/(1 + canstUo). In a similar fashion the nonlinear corrections in 3d prevent 
the divergence of 3d at the critical temperature (see the expression for the t-matrix below). 

Self- energies in t-matrix approximation. The solution to Eq. |3] at q = in the present 
ansatz may be written 
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, . . _ U$g + U T 3 5 d g + (1 - 5 d g 2 ) cos20ri 

TkM " i-uhl +5d FWMW ' ( ) 



where go and #2 are the components of the momentum integrated Green function, g = 
(1/2) J2k TrG(k, uj) and g 2 = (1/2) J2k Tr r i cos 20 G(k, a;). This solution is not exact but 
is obtained under the reasonable assumption of the smallness of higher order scattering 



terms which have no inhomogeneous driving term. ]TJ] The disorder-averaged self-energy 
is now defined in the limit of independent scattering centers to be S(k, uj) = 71^^(00), and 
determined self-consistently with the averaged G via the Dyson equation, G~ l = uj — — 
AfcTi — S(k, u) = uj — £ki~3 — AjcTi. The first term in Eq. |^ is the result obtained in the 
usual dirty d-wave theory for arbitrary scattering phase shift So = — cot _1 (l/./Vot/o)- For 
strong scattering, Sq ~ tt/2, a resonance occurs at or near the Fermi level, as can be seen by 
examining the corresponding denominator; this leads to the finite density of states at the 
Fermi level and "gapless" behavior, as has been discussed extensively in the literature. 

The denominator in the second term, due to off-diagonal scattering, leads to a similar 
resonance. To estimate the position of this feature, we examine the clean limit, in which g 2 — 
-itN [2/ii + i(uj 2 /A 2 )\og(AA /uj)} and g Q ~ -ttN [uj log(4A /cu) + iu/A ] for uj/A < I. 
The resonance in the second term of Eq. [5] then occurs when uj/Aq = [— 2/tt + (irSdNo) -1 ] = 
df. Our self-consistent determination of 3d given above shows that as Uq — > 00, Cf ~ 0.2, 
so the resonance from the off-diagonal channel never occurs at the Fermi level. It is easy 
to check that on resonance the imaginary part of the denominator is of the same order as 
the real part, up to log corrections. We therefore expect that the resonance will modify the 
low-energy behavior of the quasiparticle relaxation time at the lowest energies in the clean 
limit. 

These notions are confirmed by the full numerical evaluation of the self-energies. In Fig. 
H] we show the imaginary part of the diagonal self-energy for various values of the impurity 
scattering rate parameter T = rii/irNQ. In the clean limit the two resonances are clearly 
distinguishable, but they merge with increasing disorder. This is because the zero frequency 



feature grows only cx 7 ~ \/TA but the feature at finite frequencies oc Y (at least initially). 



Thus, above a certain disorder we no longer expect off-diagonal scattering to qualitatively 
modify low-temperature transport. There should, however, be novel temperature-dependent 
effects at small disorder. In addition, the disorder dependence of any quantity which is 
sensitive to all energy scales, as, e.g. the T = superfluid density discussed below, may be 
substantially modified. 
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FIG. 2. ImSo(w) for several values of T, with and without OP scattering. The resonance at 
finite energy is a direct consequence of the additional OP scattering. The zero frequency resonance 
is qualitatively unchanged. 



Microwave conductivity and penetration depth. In the limit of vanishing external fre- 
quency Q = 0, the quasiparticle conductivity is given by 

,2 ,oo 



o-(t, n = o) 



ne 



m 



du . 

oo \ OUO 



S(co,T), 



(6) 



where n is the electron density, e and m the electron charge and mass, respectively, and 



£ 2 



A" A' 



16 



(7) 



where £± = iyu) 2 -. — A 2 .^., the subscripts ± indicate evaluation at w ± z0 + , and real and 
imaginary parts are denoted by ' and ", respectively. The first term in Eq. [7] gives rise to the 
"universal" T — > conductivity o"oo, while the second term determines the low-temperature 
behavior. 
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FIG. 3. Low-T microwave conductivity at various disorder T, with full T-range in inset. Note 
the quasi-linear behavior of the data with OP scattering. The conductivity without OP scattering 
is always quadratic in T at low temperatures. 

In Figure |3], we show the results of a numerical evaluation of ^| with and without OP 
scattering. In the absence of OP scattering, a ~ T 2 for almost the entire region (the 
conductivity is measured in units of the "universal" zero temperature value <Jqo). The values 
of r and A Q = 3T C have been chosen by fitting the low temperature penetration depth of 
recent experiments on YBCO single crystals. |8|||. This result for the conductivity is at 
odds with the observed dominant linear behavior observed in the same samples. In contrast, 
the result including OP scattering displays quasi-linear behavior above an energy scale that 
is a fraction of 7 due to the new resonance. This becomes more obvious if we lower the 
parameter T in the case of additional OPS, so that the conductivity at higher temperatures 
is in rough agreement with the standard theory and the experimental data (not shown in 
the figure). The additional feature in the self energy can increase the conductivity by more 
than a factor of two at the appropriate temperatures, and therefore resembles much more 
the dominant linear behavior of the experimental data. 

In the inset to Fig. |3], we show the full temperature range of the conductivity using the 
model of Ref. || for the inelastic scattering. The peak in a and subsequent drop at higher 
T is due to the competition between impurity and inelastic scattering. 

The superfluid density and penetration depth are given by evaluating the imaginary part 
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of the conductivity at oo = 0, q — > if we ignore nonlocal corrections which can 
be important for T<1K in very clean samples: WB\ 



oj r dd)„ A 2 , 







p s = - / duotanh — / —Re ^ — - (8) 

'- 2T J 2tt {Co 2 - A 2 ) 3 / 2 v ; 



Because we have accounted for a new source of scattering, we expect a faster depletion of the 
superfluid density for a given impurity concentration. In Figure |], we show the dependence 
of p s at T = on disorder with and without OP scattering. The comparison is complicated 
by the fact that two cases have different critical concentrations ( the concentration at which 
T c is suppressed to zero temperature). The T c (and consequently the critical concentration) 
including OP scattering is much suppressed. This is an artifact of our treatment of the OP 
scattering which is invalid when the range of the OP suppression becomes large (close to 
T c or the critical concentration). In order to make the comparison meaningful we scale the 
result of the OP scattering case so that r c = n c J nN is the same as in the standard case 
without OP scattering. 

The increase in the initial slope of the p s suppression over the usual dirty d-wave approach 
is about 30-40%. This is less than the difference noted in Ref. @ in damaged YBCO films. 
We are not aware of measurements of the absolute scale of p s in which the samples are 
systematically disordered by, e.g. Zn substitution for Cu. Such experiments would be very 
useful to confirm the importance of order parameter scattering for bulk properties. 

Low temperature properties dependent primarily on the density of states (DOS) should 
be largely unaffected by the new source of scattering. For example, the T — > T 2 crossover 
with disorder in the penetration depth of a d-wave superconductor []T7J is controlled by the 



finite residual density of states at zero energy, N{oo — > 0), that arises due to impurities. As 
this residual DOS is determined primarily by the resonance of the self energy at the Fermi 
level in the unitarity limit, little change from the standard theory is observed. 
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FIG. 4. Superfluid density p s vs. V with & w/o OP scattering. In order to eliminate trivial 
effects due to the suppressed critical concentration in the OP scattering case we scale the corre- 
sponding data so that the T c are the same T c ~ 0.82. Note the increased initial slope in the OP 
scattering case. Close to T c our approximate treatment of OP scattering becomes invalid. 

Conclusions. The strong suppression of the order parameter around nonmagnetic impu- 
rity sites is an important qualitative difference between d-wave systems and classic s-wave 
superconductors, where such suppressions are small. The spatial structure of the quasi- 
particle states and condensate in the neighborhood of an impurity in the Cu-0 plane will 
soon be probed with STM as has been achieved for the similar vortex problem. We have 
posed the complementary question of how disorder in the off-diagonal channel induced by 
impurities modifies bulk average properties of a d-wave superconductor. By making the 
plausible assumption of local off-diagonal scattering, we have constructed a simple theory 
able to describe these effects. We have shown that the measured quasi-linear temperature 
dependence of the microwave conductivity at low temperatures can be obtained within our 
theory, as a result of a finite-energy scattering resonance which arises in the unitarity limit. 
We have also shown that the suppression of the T = superfluid density with disorder is en- 
hanced over the usual theory, although the temperature dependence is only weakly affected. 
Systematic measurements of the disorder dependence of A(T — > 0) or p s (T — > 0) will aid in 
establishing the validity of our picture. 
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